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We study the adiabatic dynamics of the Lipkin-Meshkov-Glick (LMG) close to its quantum critical 
point by linearly switching the transverse field from an initial large value to zero. We concentrate 
our attention on the residual energy after the quench in order to characterize the level of diabaticity 
of the evolution. As a function of the characteristic time of the quench r we identify three different 
regimes. For fast quenches the residual energy E Tes is almost independent on r. For slower quenches 
a second intermediate region appears in which a power-like decay emerges with E ICS ~ r~ 3//2 . By 
further slowing the quench rate, we find a third large-r regime characterized by a different power- 
law, E leB ~ r~ 2 . All these findings can be accounted for by means of an effective Landau-Zener 
approximation for the finite size LMG model. We complete our description of the adiabatic dynamics 
of the LMG model through the analysis of the entanglement entropy of the evolved state. 



PACS numbers: 



I. INTRODUCTION 



Understanding the adiabatic quantum dynamics of 
many-body systems is central to many areas of physics 
and information science. In adiabatic quantum compu- 
tationi, or in quantum annealing (see [2| for a review) 
the ultimate goal is to find the ground state of a com- 
plex system by adiabatically transforming the underlying 
Hamiltonian. Indeed any quantum algorithm can be effi- 
ciently implemented through the adiabatic evolution of a 
system from an initial exactly known state toward a un- 
known final one in which encodes the answer to the spe- 
cific computational taski Once the time scale on which 
the Hamiltonian is varied is large compared to the typical 
inverse spectral gap of the system, the quantum adiabatic 
theorem^ ensures that if the system was prepared in the 
ground state of the initial Hamiltonian, at the end of the 
evolution the quantum state will be the ground state of 
the final Hamiltonian. The bottleneck to the speed at 
which the algorithm is performed is thus given by those 
places where the instantaneous Hamiltonian has a spec- 
trum where the gap closes in the thermodynamic limit, 
i.e. as the number of qubits increases. If the minimum 
gap closes faster than a power of the number of qubits 
then the corresponding computational task is intractable. 

The closing of a gap between the ground state and 
the first excited level in the thermodynamic limit is a 
distinct feature of second order quantum phase transi- 
tion. It is responsible for the critical slowling dowr£, 
and the evolution becomes necessarily not adiabatic. 
The problem of adiabatic dynamics close to a critical 



point, and the consequent defects formation, has orig- 
inally born in the study of phase transitions in the 
early universe^ 7 -. The recent extension to the quan- 
tum case^ has stimulated an intense theoretical activ- 

^- ,,10.11,12,13,14,15,16.17,18.19,20,21,22,23.24,25,26,27,28,29 

In the search for a deeper understanding of the loss 
of adiabaticity on crossing a quantum critical point an 
important role is played by exactly solvable models. The 
study of non-equilibrium many-body system is indeed a 
formidable task and the help of a tractable exactly solv- 
able system can be of great help in testing approximate 
approaches, besides being of interest in itself. Most of 
the work done so far in this direction concentrated on 
one-dimensional quantum systems with short range in- 
teraction. In this paper we would like address a comple- 
mentary limit, i.e. a model with infinite coordination (in 
the thermodynamic limit), but still amenable of an exact 
solution: the Lipkin-Meschkov-Glick model (LMG). First 
introduced by Lipkin, Meschkov and Glick4£ in the con- 
text of nuclear physics, it was then adopted by the con- 
densed matter community as paradigm of an infinitely 
coordinated solvable system^!. The result of a sudden 
quench in this model was recently discussed ir>22, here we 
present results in the opposite regime in which the sys- 
tem is dragged adiabatically through the critical point. 
As it will be shown in the following, although the phase 
transition is of mean field nature, the dynamics leads to 
non-trivial results. 

The paper is organized as follows: In Sec. [II] we in- 
troduce the Lipkin-Meschkov-Glick model and briefly re- 
view its properties which are important for the purposes 
of this work. In the same section we also discuss how wc 
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solve numerically the dynamics, Sec. Ill Al and the observ- 
ables used to quantify the departure from the adiabatic 
ground state, Sec. Ill Bl and Sec. Ill Cl In this work we use 
the residual energy (the excess energy as compared to 
the adiabatic limit), the incomplete magnetization (the 
deficit magnetization as compared to the adiabatic limit) 
and the entanglement entropy. The numerical results to- 
gether with the corresponding scaling arguments are pre- 
sented in Sec. lIIII In the final section we present a critical 
assessment of our findings. 



II. THE MODEL 

The properties of the LMG model have been thor- 
ou ghly s crutinized in the literature (see e.g. [33ll34ll35ll36l . 
H^HHIiliOlilllillillillil and references therein). Be- 
low we briefly recall few results that are relevant to the 
present paper. The LMG Hamiltonian describes a system 
of spins (1/2 in this work) interacting through an infinite- 
range exchange coupling and immersed in a transverse 
field. Assuming that the field is directed along the z- 
direction the Hamiltonian can be written as 



TV 



H = -jf E(^ 5 J + - r E s ? 



(i) 



i<j 



where N is the number of the spins in the system, Si are 
the Pauli operators, 7 is the anisotropy parameter and T 
is the transverse field. By introducing the total spin oper- 
ator S = the Hamiltonian can be rewritten, apart 
from a additive constant, as H = — jj[S 2 + 7<->«] — FS Z . 
The Hamiltonian hence commutes with S 2 and does not 
couple states having a different parity of the number of 
spins pointing in the magnetic field direction: [H, S 2 ] = 
and [-ff, Ili <Sf] = 0- I 11 the isotropic case 7=1 also the 
z-component S is conserved, [H,S Z ] = 0. 

In the thermodynamical limit the LMG model under- 
goes a second order quantum phase transition at T c = 1 
characterized by mean-field critical exponents 3 ^. The 
magnetization in the x-direction (or in the xy-plane, for 
7 = 1) vanishes when T — > 1 _ as 



(i_ r 2)i/2 r <! 



(2) 







r > 1 



for all values of the anisotropy parameter 7. For T > r c 
and for any 7 the ground state is non degenerate; while 
for r < r c it is doubly degenerate in the thermodynam- 
ical limit for any 7 7^ 1, signaling the breaking of the Z2 
symmetry. The gap vanishes at the transition as 



A= [(r-lXr-7)] 1 / 2 forr>l. 



(3) 



For any finite N both the magnetization and the gap 
are modified (as any other physical observable). The 
finite size scaling behavior is available in literature in 



all the relevant regimes (see, e.g., f33H4Cj ). The devia- 
tion from the thermodynamic limit for the gap 5 An = 
Ajy — A and the magnetization 5m n = — m scale as 



5A N ~ N- 1 
5A N ~ 
A(N) ~ e~ aN 



5m N ~ N~ x / 2 T > 1 
5m N ~ TV" 1 / 3 T = 1 
Sm N ~ AT 1 T < 1 



(4) 



for 7 > 1 (where a is a constant) and 



5A N 
5A N 
5A N 



N- 1 
N- 1 
N" 1 



5m N 
SrriN 
5m n 



AT-i/2 r > 1 
N -i/2 T = 1 

N' 1 r < 1 



(5) 



for 7 = 1, respectively. The scaling behavior of the gap is 
important in order to distinguish the various dynamical 
regimes in the adiabatic annealing. It is however impor- 
tant to stress at this point that the equilibrium gap is 
not necessarily the one responsible for the loss of adia- 
baticity. As we will see in the following section, due to 
the parity conservation the relevant gap for the dynamics 
is different from the equilibrium one (although with the 
same scaling behavior). 



A. Adiabatic dynamics 

The adiabatic dynamics is implemented by changing 
the external transverse field from an initial value T ^> 1 
at tj n = —00, where the ground state of H(ti n ) is com- 
pletely dominated by the transverse field term with all 
the spins aligned along the +z direction, to T = 0, where 
the ground state is ordered in the xy plane. The an- 
nealing time is characterized by a time scale r. More 
specifically we consider the often in this type of 

problems, to reduce the magnetic field linearly in time 



T(t) = -t/r for te(-t in ,0] 



(G) 



with t{ n 2> T. 

The problem we want to discuss is further simplified by 
the following observation. The initial state, the ground 
state of ff(iin), belongs to the sector of maximum spin 
S = N/2. Since S is a constant of motion it is sufficient 
to restrict our attention to this subspace only. From now 
on we assume S = N/2 (for simplicity we consider 
even). In the basis \N/2,S Z ) (S z = -N/2, N/2), the 
Schrodinger evolution of the state 

N/2+1 

m))= E u 2j -i(t)\N/2,-N/2-2 + 2j), (7) 
3=1 

amounts to solving the following set of coupled equations 
.du 2] -i 

k 



dt 



E Aj\k U 2fc-l(t) 



(8) 



The odd amplitudes [JV/2, -A/2 - 1 + 2j) do not couple 
because of parity conservation. In Eq.© A is a (Af/2 + 
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1) x (N/2 + 1) symmetric matrix whose non-zero entries 
are given by 

1 

- - l) a -N/2-2+2 a-N/2+2j-l 

= + T)[ a -Af/2-3+2j + a ~N/2-2+2j] 

N 1 
-r(- y -2 + 2i) + -(l+ 7 ), (9) 

in terms of the usual angular momentum raising operator 
matrix elements: 



N ,N 
y(y + l)-iO + l) 



-,1/2 



(10) 



Special values have the boundary terms of A, given by: 
A ltl ■ 



L (l + 1 )a 2 _ N/2 -T(-^) 



1 AT/2+l, N/2+1 



4N 

+ J(l + 7) 
_ 4iV (1 + 7)a ^ 2 - 1 

-r(^) + l(i + 7 ) 



(ii) 



As discussed in Botet et a^i, the previous definition dif- 
fers from that of the spontaneous magnetization; how- 
ever, it is more amenable for finite size systems and it 
reduces to the spontaneous magnetization in the ther- 
modynamic limit. 

Since we are dealing with a model where the coupling 
has an infinite range, the incomplete magnetization is an 
appropriate way to characterize the loss of adiabaticity. 
In this case a correlation length characterizing the typi- 
cal distance between defects, along the lines followed for 
short range models, cannot be introduced. 

In the Ising limit, 7 = 0, at T(t = 0) = 0, the 
residual energy and the incomplete magnetization are re- 
lated, as they both depend only on the average value 
(4>(t = 0)\S 2 \^p(t = 0)): The residual energy per site can 
be expressed as 



Ere 



N 



.^(4 = 0)15^(4 = 0)) + ^ 



the incomplete magnetization is given by 



= 1 - 



W = 0)\S*m = 0)) 



N 2 



(15) 



(16) 



The equations ((5|) were integrated via standard numerical 
methods with initial conditions given by the amplitudes 
of the ground state of H(t = t; n ). 



B. Residual energy and incomplete magnetization 

A natural way of quantifying the degree of adiabaticity 
of the evolution is to measure the residual energy, defined 
as 



E les — Ef\ n — E. 



(12) 



where E gs is the ground state energy of H(tf\ n ), and 
E&a = (ip(tfm)\H(t fill )\'ip(t &n )) is the average energy of 
the final time-evolved state \ip(tf\ n )). Obviously Ea n , and 
hence -E7 re sj depends on the annealing time r; the slower 
the evolution the closer the final energy to E gs , hence the 
small the residual energy. 

An alternative way of quantifying the degree of adi- 
abaticity of the evolution is in terms of the incomplete 
magnetization in the final state, defined by 



- m(t) 



(13) 



where m gs is the static magnetization of the ground state 
of r = and m gs is the average magnetization of the final 
evolved state. Following Botet et the magnetization 
m has been defined through 



N 2 



(14) 



where the expectation value can be taken either the 
ground state, for m gs , or in the evolved state, for m(t). 



C. Entanglement entropy 

In addition to the previous observables, it was recently 
shown that important information of the lack of adia- 
baticity in the system can be acquired by analyzing the 
entanglement entropy 6 * 12 ' 15 . The study of entropy and 
other measures of entanglement has been recently inten- 
sively studied to characterize both equilibrium and non- 
equilibrium quantum many-body systems (see [46l | for a 
review). In the case of the LMG model the ground state 
entanglement entropy was studied i n 47 ' 48 ' 49 ' 50 ' 51 . In the 
present work we study the time evolution of S during an 
adiabatic evolution. 

Given a bipartition of the system in L and N — L spins 
the entanglement entropy associated to the reduced den- 
sity matrix of one of the subsystems, say pl, is defined 
as 



-Tr(p L log 2 p L ) 



(17) 



The entropy Sl measures the entanglement between the 
L spins and the rest of the system. 

The entanglement entropy is straightforwardly evalu- 
ated by noticing that, being the states |5 = N/2,S Z ) 
symmetric under any permutations of the sites and be- 
ing the maximum value of the total spin achievable only 
with the maximum value of the spin in each subsystem, 
the following decomposition holds^ 8 -: 



\N/2,S Z 



= x>; /2 iw-l/2) 

® \{N - L)/2,n-l- (N-L)/2) . 
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FIG. 1: (Color online) Residual energy E les (t) versus t for 
a given instance with TV = 32, 7 = of the LMG model at 
different values of r. The solid lines are the lowest-lying in- 
stantaneous spectral gaps as a function of F. The red-dashed 
line is the best fit to the lowest gap used to calculate the 
Landau-Zener transition rates. 



In the previous decomposition n and I indicate, respec- 
tively, the number of up-spins in the system and in the 
partition which defines the L sites. The coefficients ap- 
pearing are defined as pi = L\(N — L)\n\(N — n)\/(l\(L — 
l)\(n-l)\(N-L-n + l)\N\). Withthc knowledge of the 
representation of the evolved state in the basis \N/2,S Z ) 
and by the using previous decomposition, it is immedi- 
ate to trace out the TV — L spins to obtain the reduced 
density matrix p^, and calculate its entropy. 



FIG. 2: (Color online) Smallest gap and dynamical gap at the 
critical point as function of the size of the system. 
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III. RESULTS 

The results presented below were obtained by integrat- 
ing numerically Eq. ([5]). We verified that, as for the initial 
time of the evolution, it is enough to consider t- m = — 5t 
for faster sweeps (1 < r < 500) and i; n = — 2t for slower 
ones. We checked (data not reported) that our results 
do not depend on the precise value of t ln . We considered 
systems up to TV = 1024 spins and annealing times up to 
r - 10 3 - 10 4 . 

Residual energy and incomplete magnetization - In or- 
der to understand the mechanism that leads to break- 
down of adiabaticity in the LMG model it is instructive 
to start with one particular example. In Fig[T]we chose 
a system with TV = 32 spins and 7 = 0, showing the time 
evolution of the residual energy for different values of the 
annealing time r. We also plot the instantaneous accessi- 
ble gaps (thick solid lines) obtained by diagonalizing the 
Hamiltonian at any given T. As one can see, as soon 
as the system loses the adiabaticity, for fast annealing, 
it starts to ramp up in energy. The characteristic time 
scale for breaking of adiabaticity is however not given by 
the equilibrium smallest gap. As noticed in the previous 
Section, the dynamics is restricted to the subspace with 
fixed total spin S = TV/2 and can involve only states with 



FIG. 3: (Color online) Residual energy per site and incom- 
plete magnetization for the LMG model with TV = 1024 for 
different values of the anisotropy parameter 7. In all cases, for 
slow enough quenches, a power-law behavior r~ 3 ^ 2 appears. 



the same parity of S z — . Hence, the first gap relevant for 
the dynamics, that we call dynamical gap, is the energy 
difference between the ground state and the second ex- 
cited state, the smallest gap being forbidden by parity 
conservation of the total spin along the z-axis. As shown 
in Figl2l the dynamical gap exhibits the same critical be- 
havior of the excitation gap^ 3 -: both close polynomially 
in the thermodynamical limit with the same dynamical 
exponent z = 1/3, A c ~ TV _Z . This is usually accom- 
panied by a polynomial-like decay of the residual energy 
with increasing annealing time r. This is indeed the case, 
for both the residual energy and the incomplete magne- 
tization, as shown in Fig[3]and, more in detail for 7 = 0, 
in Fig[4] The behavior appears to be qualitatively inde- 
pendent on the value of the anisotropy parameter 7 for 
7 < 1, sec Fig El this was expected due to the fact that 
the minimum gap has the same large-N behavior irre- 
spective of the anisotropy. In the following only the case 
7 = will be discussed. 

Inspection of FigH] reveals three different regimes. For 
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FIG. 4: (Color online) Residual energy per spin as function 
of t for 7 = compared with different power-law behaviors. 



fast quenches the dynamics involves almost all the levels, 
see, e.g., Fig[T]for r = 5. The residual energy per site is 
near to its maximum and shows very little dependence on 
the size of the system and on the annealing time r. For 
larger values of r, a second intermediate region appears 
in which a power- like decay emerges, with E rcB ~ t~ 3 / 2 . 
Finally, by further slowing the quench rate, a third largc- 
t regime characterized by a different power-law, E Tes ~ 
t -2 , emerges. We briefly discuss the emergence of the 
last two regimes by means of a Landau-Zencr approach 
adapted to the present problem. 

The argument follows closely the one given in£. The 
probability of exciting the system into the first excited 
state, obtained by the Landau-Zener formula 



P, 



LZ 



-qA 2 t 



(18) 



with a = 7r/4, gives a lower bound to the true transition 
probability, as it ignore the transitions to all the other 
excited levels. Using the scaling of the critical point gap 
with the number of spins, A ~ iV -1 / 3 , it is possible to 
determine maximum system size for a defect-free quench 
once the probability for this to occur is fixed to the value 



P ex . This gives: 



InR. 



3/2 



free 



T 3/2 



(19) 



One can consider 1/Nf ree as an estimate of the fraction of 
the flipped spins after the quench. The residual energy 
per site in the LMG model can then be evaluated to be 



1 



N N 2 N h 



N „ T const. 
N 



r 3/2 



(20) 



This simple estimate is in good agreement with the nu- 
merical data in the intermediate regime of FigHJ 

For short range models the same power law of Eq. (j2"0|) 
can be also derived by determining the spatial scale over 
which defects occur—. We tried to apply the arguments of 
Zurek et a& to the LMG model by identifying the correla- 
tion length with the coherence number introduced in [3l| . 
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FIG. 5: (Color online) Comparison between the excitation 
probabilities as function of F of the LMG model with N — 32 
and of its effective LZ approximation for different values of r. 



The procedure we followed, however, does not lead to the 
correct exponent. We have reasonable confidence that 
the failure in obtaining the correct scaling with this sec- 
ond method may be related to the above identification 
and the consequent definition of defect density. It would 
be interesting to find the correct argument in order to 
extend the approach by Zurek et a& or Polkovniko^ to 
infinite range models. 

Effective two-level approximation - As already men- 
tioned before there is, for slower quenches, a further 
crossover to a different power-law. Can one explain also 
this behavior by using a Landau-Zencr argument? To 
this end, it is important to refine this comparison and to 
understand to which extent the dynamics of a many-body 
system described by the LMG model can be described 
by two (many-body) levels. In general, in a many-body 
system there will be a number of avoided crossings and 
multiple LZ transitions, including interference between 
them. Only when a single avoided crossing is dominant 
and well separated from the others a two-level approxi- 
mation is appropriate. A detailed analysis of this issue 
is summarized in Figs[5] and [S] where we show the case of 
N = 32 as an example. Our analysis starts by extracting 
the best dynamical minimum gap and adapting to it the 
following two-level Hamiltonian: 



H 



LZ 



-n LZ (T-T ) A LZ 

A LZ n LZ (r-r ) 



(21) 



In the effective Landau-Zencr problem ^l^z, l\lz and 
To are the fitting parameters and T = —t/r. In Fig[T] 
the dashed line represents the instantaneous gap of the 
Hamiltonian (|2"Tj) suited to the case N — 32. From here 
we compare the results of the full LMG model with those 
obtained using LZ theory. As shown in FigO the exci- 
tation probability in the LMG model for slow enough 
quenches coincides with that of the effective LZ prob- 
lem. It appears that this approximation is good also in 
the estimate of the asymptotic value of the probability 
for 10 < r < 100. Deviations come predominantly from 
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FIG. 6: (Color online) Excitation probability as function of 
r for LMG systems (circles) of two different sizes compared 
to that one of the respective LZ-effective models (line) for 
different final times. For the LZ models the first two terms 
of the Vitanov approximation have been used. For the case 
LMG N=32, the probability of exciting the first level is also 
presented (triangles). 



the more enhanced oscillations of the post crossing re- 
gion in the LMG model. For larges r's the asymptotic 
value obtained from the effective two-level system gives 
a very poor approximation to the actual data. This can 
be traced back to the presence of further crossings which 
are obviously neglected in the two-level approximation. 
In our LZ scheme this can be effectively corrected by 
approximating the LZ crossing probability to the time 
before the next level crossing comes into play. This is 
explained below 

As found by Vitano\s£, it is possible to define the du- 
ration of a single LZ as the time required by the probabil- 
ity for jumping from zero to its asymptotic value, linearly 
and with the slope calculated at the crossing point. Using 
r as time-scale one can write: 



P(oo) 

P (Tcross) 



(22) 



This time turns out to be exponentially divergent with 
t for large r— . This means that for slow quenches 
consecutive LZ transitions are not independent. In a 
first crude approximation, we can guess that the con- 
sequence of this is simply to stop the probability from 
relaxing towards the asymptotic value when the sys- 
tem has reached the second crossing. The presence of 
a power-law regime ~ t~ 2 for extremely slow dynamics 
is a clear consequence of the finite duration of the evo- 
lution. In the original works by Landau and by Zencr, 
the final time is supposed to be tf = oo; here the evolu- 
tion is stopped at Tf = — tf/r = for the LMG model, 
and at t(LZ)f = —T{LZ)f • T for the effective LZ, with 
= — r . An accurate analysis of the finite-time 
Landau-Zener model (FTLZ) has been done in Ref. [53| . 
where it is shown that the transition probability reads, 



P 



(FTLZ) 



(r) ~ Plz(t) 



(1 - 2P lz (t)) 



1 fiA 4 -X, 



-(i 



\3 
(LZ)f) 



1 lz y-2 

' 



(23) 



with Plz{t) = e~ 7 ' A ~ LzT / nLZ . As it can be immediately 
seen from the previous equation, by sending the final time 
to infinity the usual LZ probability is recovered. The 
crossover rate f to the r -2 scaling is obtained by equating 
the two terms on the r.h.s. of Eq. ([23| . In the limit -(1 + 



r? izu ) 3 / 2 » 1 the crossover time is approximated 



ALT (£2)/ 
by 



LZ 



4A 2 /I , ®?T7 

^^lz (1 H ^ 



F 2 W 2 

lz (LZ)f> 



(24) 



In Fig[6]we compared the excitation probabilities of LMG 
systems of different sizes with their singlc-LZ approxi- 
mations. The probabilities for the effective models are 
evaluated for three different final time: Tf = 0,Fi,r2, 
where the last two are the positions, respectively, of the 
minimum gap between the ground state and the second 
excited level, and the minimum gap between the first and 
the second excited levels. As it can be seen, the agree- 
ment is quite good and one can reproduce in this way 
also the regime with the r -2 behavior. 

Entanglement entropy - We finally would like to discuss 
the behavior of the entanglement entropy which was al- 
ready used as a tool to characterize adiabatic many-body 
dynamics in Refs. (l2lfl5l | . Our results are summarized in 
Fig|7] In the left lower panel the entanglement entropy, 
for a block of size L = N/2, of the state evolved down 
to Tf — is plotted as function of the quench time r. 
For fast quenches, r — > 0,the state does not evolve (it re- 
mains in a nearly factorized state), thus the entanglement 
necessarily tends towards zero. For very slow dynamics 
t — ► oo, since we are dealing with finite systems, the 
evolution eventually becomes adiabatic and the entan- 
glement picks up the value it assumes in the final ground 
state, S gs (T f = 0) = 1, independently on the subsystem 
size^. Between this two limiting behaviors, the entropy 
reaches a size-dependent maximum at an intermediate 
value of r. An interesting feature is that the presence 
of a finite minimum gap can be easily connected with a 
time scale for the decaying of the entanglement. A pos- 
sible choice for this time scale consists in selecting the r* 
at which the entropy has reduced by half the value of the 
its peak respect to the slow quench limit 



Sn/2( t *) 



(S n 
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1 
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In the upper panel of Fig. ([7]), r* determined in this way 
is shown as a function of the system size N. For large TV, 
a power-like behavior emerges with an exponent ~ 0.66, 
hinting at a relation 
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FIG. 7: (Color online) Left lower panel: entanglement entropy 
of a block of L — N/2 spins as function of the quench time 
t. Right lower panel: entanglement entropy of a block of 
L — N/2 divided by its maximum value as function of the 
rescaled variable t/N 0,66 . Upper panel: the time scale t* , 
see the text for the definition, as function of system size N. 



In the lower right panel of Fig. ([7]), the entanglement 
Sjv/2 divided by its maximum value S mayi is plotted as 
a function of the rescaled variable t/N 0S6 , showing, for 
large systems (N > 128), a collapse of all data on the 
same curve. Note that Eq. (|2l)|) expresses exactly the 
same energy-time relation found in the usual LZ system, 
see Eq. tfTB)) . so that the correspondence stated in previous 
sections is again supported. 



evaluate the extent of diabaticity of the evolution. The 
dynamics is restricted to a subspace of definite total spin 
and parity of its projection along the z-axis, due to the 
symmetries of the Hamiltonian. Results appeared to 
be qualitatively independent of the value of the XY- 
anisotropy parameter 7, except for the fully isotropic 
XX case at 7 = 1, where the further conservation of S z 
plays an important role. By starting the evolution in the 
ground state for very large values of the transverse field 
r, and then reducing T(t) linearly to zero, three regimes 
in the residual energy are identifiable: the first one, cor- 
responding to fast quenches, is strongly not adiabatic, 
involves transitions from the ground state towards many 
excited states and is characterized by a residual energy 
near to its saturation value. In the intermediate regime, 
the lowest critical dynamically accessible gap starts dom- 
inating the evolution, inducing a residual energy per site 
that decays in a power- like manner, like r -3 / 2 . The third 
large-T region, where the residual energy decays like t~ 2 , 
is understood by taking into account the presence of ad- 
ditional level crossings. In the effective Landau-Zcncr 
description we used in this paper this results in the re- 
quirement to consider a finite-time Landau-Zener sweep. 
As show by Vitanov a finite-time sweep leads to a polyno- 
mial (in t) contribution to the LZ transition probability 
which is dominant for very slow sweeping rates. Notice 
that this t -2 regime, usually described as the general de- 
viation from adiabaticity deriving by the adiabatic the- 
orem for very slow evolutions^, emerges here in an al- 
ternative way through the parallelism with an effective 
FTLZ model. 



IV. CONCLUSIONS 

In this paper we have studied the adiabatic quantum 
dynamics of the LMG model in a transverse field across 
its quantum critical point. We focused our attention on 
the residual energy after the quench analyzing its be- 
havior as a function of the annealing time, in order to 
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